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Abstract: We study cosmological perturbations produced by the most general two- 
derivative actions involving two scalar fields, coupled to Einstein gravity, with an arbitrary 
field space metric, that admit scaling solutions. For contracting universes, we show that scale- 
invariant adiabatic perturbations can be produced continuously as modes leave the horizon 
for any equation of state parameter w > 0. The corresponding background solutions are un- 
stable, which we argue is a universal feature of contracting models that yield scale-invariant 
spectra. For expanding universes, we find that nearly scale-invariant adiabatic perturbation 
spectra can only be produced for w « — 1, and that the corresponding scaling solutions are 
attractors. The presence of a nontrivial metric on field space is a crucial ingredient in our 
results. 
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1. Introduction 

Single-field, slow-roll inflation has become the dominant paradigm for generating primordial 
density perturbations in the early universe. Despite its successes, it is still important to 
search for alternative models, if only to determine whether the predictions of inflation are 
truly unique. Any proposed alternative must explain the nearly scale-invariant adiabatic 
primordial density fluctuations that seed large scale structure formation, and whose imprint 
is seen directly in the cosmic microwave background. 

Models with a contracting phase, where density perturbations are generated before the 
traditional big bang, provide an interesting class of alternatives to inflation. Proposals of 
this type include the recent ekyprotic and cyclic models [1, 2], and the older pre-big bang 
scenario [3]. The ekpyrotic and cyclic models create primordial cosmological perturbations 
during a w > 1 contracting phase, which has some similarity to an expanding inflationary 
phase [4, 5, 6]. However, it is well known that a single scalar field in a contracting w > 1 
phase cannot produce nearly scale-invariant curvature perturbations, for ( is not sensitive to 
the nearly scale-invariant growing mode of the Newtonian potential, but only to the decaying 
mode [7, 8, 9, 10, 11]. Several approaches to dealing with this problem have been proposed. 
One can analyze the perturbations in higher dimensions [12, 13, 14] in the hope that higher- 
dimensional effects will mix the growing and decaying modes. Alternatively one may hope 
that during the evolution through the bounce from contraction to expansion a similar effect 
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occurs [15, 16, 17, 18, 19, 20]. Or, one can posit that a spectrum of scale-invariant isocurvature 
modes are subsequently converted into adiabatic modes [21, 22, 23, 24, 25], an idea which has 
also been applied in expanding universes [26]. In still other models, fluctuations are adiabatic 
and scale-invariant at horizon crossing [27]. 

In this work we describe a mechanism that generates scale-invariant spectra in both 
expanding and contracting universes. We consider the entire family of scaling solutions with 
two scalar fields, including a nontrivial metric on field space. These solutions serve as models 
of realistic systems in which the equation of state parameter w is slowly varying, and are 
useful tools for understanding the predictions of inflationary- and ekpyrotic-type models. We 
find that using our new two-field solutions it is possible to construct contracting universe 
models in which a scale-invariant spectrum of adiabatic fluctuations is continuously produced 
at horizon crossing for any equation of state parameter w > 0. The price is an instability 
in the background solution, which we show is universally present in all two-field contracting 
universe models that produce scale-invariant spectra. There are also inflationary solutions 
with nearly scale- invariant adiabatic perturbation spectra, but only when w « — 1. The 
inflationary background solutions are all stable to small perturbations. Within the class of 
two scalar field models that admit scaling solutions, both the stability of the background 
solution and the spectral indices of long-wavelength perturbations are determined by only 
two parameters. The first is the equation of state parameter w for the background solution, 
which we here parameterize by c, with c = y/3(l + w). The second parameter, which we call 
A, is constructed from the first and second derivatives of various functions that appear in the 
general action for scaling solutions. 

We will justify our claims with detailed calculations presented in the following sections. 
Nonetheless, the essence of our results can be understood using simple arguments. The first 
important fact is that the behavior of fluctuations on sub-horizon scales, with comoving 
wavenumbers k > aH, is determined by the behavior of a(rj) alone. This is because the 
action for the curvature perturbation £ is of the form (here rj is conformal time and ' denotes 
conformal time derivatives) 



= j C ^ft{iC k ?-k\l + ---) dr,, (1.1) 



where we have specialized to a single Fourier mode Ck(v) °f comoving wavenumber k. Con- 
sequently, the standard boundary conditions on Qkijf) are such that when well within the 
horizon each (positive-frequency) mode is 

-ikr; 

^" ) = 77S^ f ° r |l " l>>L (L2) 

For scaling solutions, the equation of state parameter w is constant, so a(rj) oc |ry| 2 /( 1 + 3tu ). At 
horizon crossing, where \krj\ = 1, the amplitude of each mode is thus 

C k (v)=Ck B with B = ^ ~ 3 ™ } when \k V \ = 1, (1.3) 
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with C a fc-independent constant. The modes will only be exactly scale-invariant as they 
cross the horizon, with B = —3/2, when w = — 1. 

In the absence of entropy perturbations, such as in single-inflaton models, it is well known 
that C' k = outside the horizon. In this case, the only way to generate a scale-invariant 
curvature spectrum as modes leave the horizon in an expanding phase is with w ~ —1, i.e. 
inflation. In a contracting universe with w > 1, £ is conserved in single-field models, so 
matching across the horizon cannot give a scale- invariant spectrum. When —1/3 < w < 1, 
one finds that ( grows outside the horizon as Cfc ~ |r/| 3(u '~ 1 )/( 1+3w ), which can be used to give 
Cfe a scale-invariant spectrum if w = [28, 29]. However, the w = background possesses an 
instability, and this mechanism is also problematic since a w = component in a contracting 
universe will be quickly overtaken by radiation and anisotropy as the scale factor shrinks to 
zero [4]. 

When entropy perturbations are present, such as in multi-inflaton models, then (' k need 
not vanish. It is then possible to continuously generate scale-invariant adiabatic perturbations 
at horizon crossing if £fc evolves in the correct way on super-horizon scales. Specifically, 
suppose Cfc evolves as 



for some exponent p. Then using (1.3) one finds that Ck will develop an exactly scale- invariant 
spectrum provided that 



If we are free to adjust p at will, the only remaining requirement is that modes are exiting 
the horizon during the time when primordial perturbations are generated. This requires that 
w > —1/3 for contracting solutions, and w < —1/3 for expanding solutions. We will see 
below that once w is fixed in the two-field scaling solutions, varying A provides the freedom 
to adjust p as required for the contracting solutions, provided that in addition that w > 0. 
A special case of this mechanism was considered in [27]. In the expanding case, the growing 
mode always has £ constant and this mechanism cannot be applied. 

The full analysis shows that even for the contacting case this mechanism cannot be 
implemented for scaling solutions when the metric on field space is Euclidean. That is, if 
there exists a field redefinition that takes the action for the scalar fields <I> a into the form 



then we will show that C is constant on super-horizon scales, despite the presence of entropy 
perturbations. The underlying reason is that adiabatic and entropy perturbations decouple for 
for scaling solutions with Euclidean field space metrics, and so the usual arguments regarding 
single-field adiabatic perturbations will apply [7, 11]. Therefore, in these cases the only 
possibility for generating scale-invariant adiabatic perturbations in is an expanding nearly 
de Sitter solution with w ~ — 1 or a contracting solution with w ~ 0. The mechanisms 
proposed recently by [22], and also employed in [23, 24], rely on mixing of isocurvature into 



Ckiv) ~ \V\ P for N| < 1, 



(1.4) 




(1.6) 
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adiabatic modes by turning a corner in field space, during which time the solution is far from 
scaling for a brief period. Here, the non-trivial metric on field space provides the crucial 
coupling between adiabatic and entropy perturbations on super-horizon scales. 

Ideally, the background solutions used to obtain these perturbation spectra would be 
stable, and so fine-tuned initial conditions would not be required. For contracting universes, 
the heuristic arguments given above hint that this will not be possible. For these solutions, 
which have w > 1, the expression (1.5) implies that we require p < 0, and so Oc grows on 
super-horizon scales. The growth of long-wavelength modes suggests an instability in the 
homogeneous degrees of freedom of the system. We will explicitly verify this instability in the 
scaling solutions by studying the fixed points of the autonomous dynamical system associated 
with the homogeneous field and metric modes. We will also argue that this instability must 
always be present if a scale-invariant spectrum is to be created in a contracting universe - 
even if it is not initially created in the adiabatic perturbation modes. Whether this instability 
is really a problem depends on the specific model under consideration. The instability would 
only be present during the time in which the background approximately corresponds to a 
scaling solution, and would in any event be naturally cut off by the bounce itself. In the 
expanding case, where —1 < w < —1/3, the expression (1.5) implies that p > and so £fc is 
either constant or decays on super-horizon scales in order to produce a nearly scale-invariant 
spectrum. This suggests that the inflationary background solutions are stable, which we 
verify by studying the associated autonomous system. 

Inflation models with multiple fields have been a fertile research topic for some time: see 
for example [30]-[65]. Some basic ingredients of the mechanism we describe here, such as the 
generation of isocurvature perturbations or non-flat metrics on field space, have been studied 
before in the context of slow-roll inflation. Our novel contribution is to discard the slow-roll 
assumption by concentrating on scaling solutions, which allows for an exact calculation of 
the spectral indices in terms of (c, A). Our results indicate that it is still necessary to be 
close to slow-roll, i.e. have w ~ — 1 to achieve an approximately scale-invariant spectrum. In 
collapsing models there is no natural slow-roll condition, and models with multiple fields that 
produce scale-invariant isocurvature fluctuations at horizon crossing, which are later converted 
to adiabatic fluctuations, have attracted a great deal of attention recently [22, 23, 24, 25]. 
In the present work we directly produce adiabatic fluctuations at horizon crossing, thanks to 
the non-Euclidean metric on field space. 

The scaling solutions that we employ in this analysis are the most general for two- 
derivative actions with two fields. We include an arbitrary metric on field space, but otherwise 
the kinetic terms are canonical. Once the equation of state parameter w is fixed for the back- 
ground solution, the relevant action is completely specified by three additional free parameters 
fi, $2 and hi-, which arise from the expansion to quadratic order of some free functions appear- 
ing in the scaling action. A remarkable feature of these scaling solutions is that their stability 
properties and the spectral index of the perturbations depend entirely on c = \/3(l + w) 
and on a single additional parameter, which we call A, that is itself constructed from /i,/2 
and /i2- This simple dependence makes it quite straightforward to study the phenomenology 
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of the two field model, as well as to include a variety of other models studied by previous 
researchers within the framework described here. In these senses we have characterized the 
natural two-field two-derivative generalizations of the single-field scaling solutions which have 
been studied in both expanding and collapsing universes [67, 5, 68]. 

In the present work we do not attempt to embed our mechanism within a complete 
model. Instead, we focus on questions of principle, and endeavor to show that it is possible to 
construct models with a simple field content and Lagrangian that will produce certain types of 
primordial perturbation spectra. We find that the combination of including a curved metric on 
field is enough to drastically modify the predictions from the single-field case for contracting 
models, but not for expanding ones. Questions such as whether the initial conditions for 
the contracting solutions are truly fine-tuned will depend on the details of specific models. 
The relative amplitudes of the curvature and isocurvature perturbations after the universe 
enters radiation domination will also depend on the details of reheating after inflation (or its 
analogue in contracting universe models). We will not consider the details of the process by 
which inflation ends in expanding models, or the big crunch/big bang transition in contracting 
models. 

As in other contracting universe models, the mechanism described here produces its own 
distinctive observational signature in the tensor spectrum. While the spectral index for scalar 
perturbations depends on the detailed coupling between adiabatic and isocurvature modes 
on super-horizon scales, the tensor spectrum does not. The tensor modes decouple from 
everything but the background metric, and so they directly probe the expansion history (and 
therefore w) during the time in which the perturbations are generated. The kinetic part of 
the action for the tensor perturbation hk is given by (1.1) under the replacement cC,k hk- 
The tensor mode hk evolves like (1.2) while sub-horizon and is constant when super-horizon, 
so by (1.3) the tensor spectral index ut is 



In contracting models, the tensor spectrum remains quite blue, with 2 > nx > 3. The more 
general solutions we study are therefore very relevant for searches for primordial gravitational 
waves in the cosmic microwave background or with laser interferometers. 

This paper is organized as follows. In Section 2 we describe the class of two scalar 
field models that admit scaling solutions. Any such model can be completely described by 
four parameters c, /i,/2, and hi- We also give explicit examples, one with a Euclidean field 
space metric, and one in which field space is curved. In Section 3 we set up and solve the 
cosmological perturbation equations for the scaling backgrounds. We show that the spectral 
indices of isocurvature and adiabatic modes depend only on c and one additional parameter 
A, which is constructed from c, /i,/2, and hi- We also show how super-horizon isocurvature 
and adiabatic modes decouple when the field space metric is Euclidean. Section 4 is devoted 
to an analysis of the stability of the background scaling solutions, which is also completely 
controlled by c and A. We conclude in Section 5. 




(1.7) 
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2. Scaling solutions 



A scaling solution is a solution to the equations of motion for which all contributions to 
the total energy density scale identically with time, keeping their fractional contributions 
constant. Although the system is evolving in time, there is a sense in which it is time 
translation invariant. To be precise, scaling solutions arise when there exists a timelike 
homothetic Killing vector 1 on spacetime, i.e. a scale invariance. 

Consider the action for a set of scalar fields <J> a with a nontrivial field space metric G a b(&) 
coupled to gravity 2 

S = j (^R-\G ah g^d^ a d^ b -vW) V=9d 4 x. (2.1) 

A scaling solution exists when there is a continuous transformation with parameter A such 
that 

^ = m 9»u^e x g^ S^e x S. (2.2) 

Any such transformation preserves the form of the equations of motion 3 . These conditions 
imply that £ a is a Killing vector on field space, and that the potential transforms as 

V -► e~ x V. (2.3) 

If we choose coordinates on field space so that the Killing direction is cp, and a 1 , . . . ,a n are 
the other directions, then we have sufficient coordinate freedom to put the metric in the form 

G ab d$ a d<$> b = /(a 1 , . . . ,a n )(dcP) 2 + 2A i (a 1 , . . . ,a n )da i dc/) + h^a 1 , . . ./JSuW, (2.4) 

and the potential in the form 

V = Vbexp(-c0)/i((T 1 ,...,cj n ). (2.5) 

Here Ai is only defined up to Ai — > Ai + J^j-, corresponding to redefinitions <fi — ► <fi + 
x((7 1 . . . a" - ), with an associated redefinition of hij. In the two field case (n = 1) we can use 
this freedom to set A\ = 0. The metric then satisfies the conditions (2.2) with a Killing vector 
given by £ a = (1/c, 0, . . . , 0). In these coordinates, the scaling solution has only <f> changing 
with time, and cr 1 ,...cr n are constants which can be taken to be zero by an appropriate 
coordinate transformation. It is then convenient to rescale 4> so that /(0, ...,0) = 1 and 
rescale Vq so that h(0, . . . , 0) = 1. The net result is that the scaling solution is just the usual 
single field scaling solution associated with the action 

Segmv = j (^R ~ \ (d</>) 2 ~ V exp(- C 0)) ^d 4 X. (2.6) 

Nevertheless perturbations around this solution will be very different, since fluctuations in 
(p and the a % will probe the terms in the Taylor expansions of f{a l ) and h(a l ) about the 
background solution. 

X A homothetic Killing vector £ a satisfies V a £f, + Vft£ a = Cg a b where C is a constant. 
2 We shall work in units where M^ 2 — 8ttG = 1. 

3 A slightly different approach to multi-field scaling solutions has been considered in [66]. 
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2.1 The general case 

The action for any two-field system that admits scaling solutions can be put in the form 

S = J Qfl-^/(<r)W) 2 ~\{daf -y exp(-c^(<7)) V=gd 4 x. (2.7) 

There is sufficient coordinate freedom to set a = in the background solution, so for linearized 
perturbations we will only need the second order Taylor expansions of the two functions 

f(a) = l + ha + f 2 a 2 + ..., (2.8) 
h(a) = 1 + K\a + h 2 a 2 + . . . . (2.9) 

This is a useful simplification since the result will depend on only a handful of coefficients. 
In fact, for a = to be a consistent solution we must have 

% = -c 2 /i/(c 2 -6). (2.10) 

This is not a physical restriction on the choice of potential and kinetic term, but instead 
reflects our choice of coordinates on field space. The coefficient Vo is degenerate with shifts 
<j) — ► 4>+C, but we will leave this symmetry unfixed. The net result is that the free parameters 
in the Lagrangian are (c, /i, / 2 , h 2 ). 

2.2 A Euclidean example 

To illustrate the above analysis it is useful to look at an example recently considered by 
Lehners et al. [22] in the context of the ekpyrotic model (see also [62, 63, 64] in the context 
of assisted inflation), of two minimally coupled fields with exponential potentials 

S = J Qi?-^(^i) 2 -i(^ 2 ) 2 -y exp(-ci$i)-Vbexp(-c 2 $ 2 )^ V=gd 4 x, (2.11) 

where we have used one of the shift symmetries <J>j — > <&j + Cj to equalize the coeffients in 
front of the exponentials. In these coordinates the scaling solution has 

<S>j(ri) = Ajlnin) + Bj, (2.12) 
with c\A\ = C2^4 2 - We can perform a rotation on field space to define 

' + (2.13) 

+ <t , (2.14) 



V4 + 4 

Cl$l - c 2 $ 2 



V4 + 4 



so that the action becomes 



S = J Qi? - ^) 2 - \(daf - Vo exp(-c0)/i(<7)) y^d 4 *, (2.15) 
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where 

? = ? + (2 ' 16) 
and we used the second shift symmetry to set 

_ 21n(c 2 /ci) 

^ " c 2 + c 2 (2 ' 17) 

which gives 



and 



V (2-18) 



/i(<j) = l + -a 2 + --- (2.19) 

so that /i 2 = c 2 /2 but /i = / 2 = ^i = 0. The scaling solution considered in [22] corresponds 
to 

^ = <Po HV) (2-20) 

and (7 = 0. This case is very special since the field space metric is Euclidean and so there are 
no metric couplings. 

2.3 A curved example 

As an example of a two-field model with a nontrivial field space metric which still exhibits 
scaling solutions, consider the following action 

( , 21) 

The field space metric is a two dimensional hyperboloid, with constant negative curvature. 
This kinetic term arises often in string compactifications, where it is more common to define 
z = tp + ib so that the kinetic term is 

(db) 2 + <W? dzdz 

A 2 6 2 [Alm(z)] 2 ' 1 " ) 

Through the redefinition b = A^ 1 e~ Ax the kinetic term is cast into another useful form 

{db) 2 + (di>) 2 



A 2 6 2 



(Sx) 2 + e 2Ax W) 2 j (2 _ 23) 



which has been studied a number of times in the past. Returning to the original variables 
(b,tp), we seek scaling solutions with equation of state parameter w, having 

b M = m = (v_\ 3A{w - 1)/il+3w) (224) 

bo Y>o VW 
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and a(rj) oc \rj\ 2 ^ 1+3w K These solutions exist provided that 



A 



a 



4A 2 1 + w 
3 (l-w)'- 

21 + t»l- 



1 + 



[V>o/&o] 2 ^ 



Al-w l + [V>oA) 



12 ' 



/3 = 



21 + w 
Al 



w 



fo_ 
bo 



l + A 



(2.25) 
(2.26) 
(2.27) 



and so if (a, (3) are freely adjustable then the scaling solutions are specified by (w, A, ipo/bo). 
The parameter ipo/bo is essentially the ratio of the kinetic energies of the two fields. With 
the field redefinition 



b = 



1> 



V> (l + C 2 e 2Ao 
i/. (l + CV Aff ) ' 



1) 



the action becomes 
S = 

where 



d 4 Xy 



m = 



h(a) 



( e -A^ + c 2 e ACT ) 2 
(1 + C 2 ) 2 : 



(1 + C 2 )e A(T 
1 + C 2 e 2ACT 



C 2 e 2Xa + 1 C 2_ 1 



(2.28) 
(2.29) 

(2.30) 

(2.31) 
(2.32) 

(2.33) 



To obtain these expressions we have applied the requirement that the scaling solution corre- 
sponds to a = and that /(0) = h(0) = 1, which fixes 

. 2 - 4C2A2 (2 34) 

As we will see below, the quantity that determines whether isocurvature and adiabatic per- 
turbations are coupled on super-horizon scales is f±, which in this case is given by 

r c 2 -r 



/i = 2A 



(2.35) 



c 2 + i. 

The variable C parameterizes the family of scaling solutions with fixed w, A. In terms of the 
ratio of kinetic energies, we have the relation 

1 



b 2 



C 



(2.36) 



For the special value C = 1, the coupling parameter f\ vanishes. This special value of C 
corresponds to ipo/bo = 0, which is the case in which there is no kinetic energy in the tp field. 
The quantities fn h\ and hi can be readily calculated from the expressions for /(<r) and h(a). 
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3. Two-field perturbations 



In this section we analyze the perturbations produced by the two field system in uniform-field 
gauge 5(f) = 0. To connect quantities in this gauge to gauge-invariant variables, we define the 
conventional gauge-invariant adiabatic perturbation variable Q (also called TZ) as 

where ip is the spatial metric perturbation, and so in uniform-field gauge £ = tjj. The fluc- 
tuation 5a is automatically gauge invariant since a = in the background. To stress the 
connection to gauge-invariant quantities we express the perturbations directly in terms of ( 
and 5a. Had we chosen to use coordinates on field space where both fields are varying (e.g. 
the $i and <E>2 considered above), a more complex projection would be required to separate 
the adiabatic and isocurvature modes (see for example [55] for an exposition of this tech- 
nique). Our analysis is valid for both contracting (ekpyrotic) and expanding (inflationary) 
solutions through an appropriate choice of the equation of state w. 

3.1 Background scaling solution 

The background equations of motion are 

3 («) 2 = 2 0,2+ r ,2+aVoe ~ c ^' (3 - 2) 

= 0" + 2 -<f>' + r 1 ¥-</>' a 1 - ca 2 V e- c ^r l h, (3.3) 
a da 



a 2 da da 
The background scaling solution expressed in terms of conformal time is 



(3.4) 



a(v) = (VM P , (3.5) 
2c 

a = 0, (3.7) 



4>0 + ^2 n^V/Vo), (3.6) 



where p = 2/(1 + 3w) and w = c 2 /3 — 1 is the equation of state parameter. The parameter 
00 is given by 

0O= c ln V 2(c 2 -6) )■ (3 ' 8) 

Note that Vq > for c 2 < 6 and Vq < for c 2 > 6. Consistency of the solution a = with 
(3.4) enforces h x = -c 2 /i/(c 2 - 6). 
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3.2 Uniform-field gauge perturbations 



The calculation of the two field perturbations is straightforward to perform in uniform-field 
gauge. The scalar metric perturbations are defined via 

ds 2 = a{7]) 2 (-(1 + 2$)d?? 2 + 2V i adx i dr j + (1 - 2C)df 2 ) , (3.9) 

where we have used the residual gauge freedom to set anisotropic parts of the 3-metric to 
zero. <1> and a are essentially the lapse and shift perturbations, respectively, in the ADM 
decomposition. They can be uniquely determined in momentum space in terms of the physical 
degrees of freedom by solving the two constraint equations 

ewj AT o 2c 2 A^ fc - 2(c 2 - 6)$ fc + (c 2 - 2)((c 2 - 2)kW(k + GyQ mm 

^0 = ^0^ $fe = _l (c 2_ 2K , (3n) 

The remaining Einstein equations are redundant as a consequence of the Bianchi identities, 
and so the equations of motion for the two physical degrees of freedom can be obtained 
directly from the two perturbed equations for the scalar fields 

,n | 4 2 2(c 2 -6)h 2/i 

Cfe + (^2^ Cfc + ^ Cfc = (c 2 -2)V " " W^Wl *' ( } 

c „ , 4 ; 2 (24/i 2 -4c 2 (/ 2 + /i 2 )) x 2c 2 /! 

fe (c 2 -2)?? ft (c 2 -2) 2 7? 2 (c 2 - 2)77 

These equations of motion arise from the quadratic action <S , (2)> obtained by expanding the 
original action (2.7) to second order and re-expressing it in terms of our perturbation variables 
Cfc, Sa k , which gives 



= / ( l(Q 2 + ^(^) 2 - - ^kV k 

(24/> 2 - 4c 2 (/ 2 + /i 2 )K 2 2/i 



c 2 (c 2 _ 2)2r? 2 6a k + (c 2_ 2)r/ ^Cfc J c a(r?) dry. (3.14) 

The nature of the solutions depends crucially on the kinetic coupling parameter f±. If f± = 
then the adiabatic and isocurvature modes decouple: this is the case when the field space 
metric is Euclidean. If f\ / then the adiabatic and isocurvature modes are inextricably 
coupled. That is, we have already used all of our reparameterizations to put the action in the 
form (3.14), so the coupling term cannot be removed. Because of the significant differences 
between the /1 = and f\ / cases, in what follows we shall consider the two solutions 
separately. 
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The inner product of two sets of modes can be inferred from the action (3.14), and is 
given in momentum space by the analogue of the Klein-Gordon norm 

<(Cfc, Sa k ), (( k ,6a k )) = ia 2 ^) ( c 2 Q( k - c 2 (>* k t k + Sa* k 5a' k 

which is defined so that a positive frequency mode has unit norm. 
Decoupled case: f\ = 

In the decoupled case we can find exact solutions for the mode functions in terms of Hankel 
functions. The two independent positive frequency modes for each field are given by 

^ = ^m^-^ H ^-^ =a (3 - 17) 



ii 



where 



and 



c 2 



v= w=ry (3 - 18) 



^( c 2_ 6) 2 + 16c 2 (/2 + (1 _ 6/c 2 )/l2) 



(3.19) 



2(c 2 - 2) 

Since the isocurvature and adiabatic modes are decoupled, they will have different long- 
wavelength spectral indices. We define the power spectra in the usual way 

Pdk,Vf) = ^(CUvf)), (3-20) 

P *(k,Vf) = ^(So&Vf)), (3-21) 

where rff is the time at which the generation process ends, such as at reheating in the infla- 
tionary case. Since £ is a constant outside the horizon, does not depend on 77/ and we find 
an adiabatic spectral index 

n s = 4-2H=4-j^|[. (3.22) 



So as for a single field the only two case which give scale invariance are c 2 = i.e. w = —I 
expanding, or c 2 = 3, w = contracting. By contrast a will in general be evolving outside 
the horizon and so its amplitude will depend on when the generation process ends. The 
isocurvature spectral index of the modes which have already left the horizon are given by 
n s = 4 — 2|z/|, which can be freely adjusted by changing the coupling constants. This fact has 
already been used in many scenarios [26, 21, 22, 23, 24] to generate scale invariant isocurvature 
modes which are subsequently converted to adiabatic ones by a second mechanism. 
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Coupled case: f\ / 

When their equations of motion are coupled, the isocurvature and adiabatic modes must 
have essentially the same dependence on rj. This implies that both modes have the same 
spectral index at long- wavelengths. We compute the spectral index by taking the usual WKB 
approximation within the horizon, matching this solution at horizon crossing with a power law 
super-horizon mode. Because our solutions are scaling this will reproduce the exact spectral 
index, but only give the amplitude approximately. 

For k 3> aH the positive frequency modes are given by their WKB form 

"i = «tfsW'' (3 - 24) 

where 

*± = JJ^y (3-25) 

B (c 2 -2)s , , 

- = - { - ' = T*c. (3.26) 

Again we have two independent pairs of solutions and we should distinguish between the two 
which are orthogonal with respect to (3.15). In fact the above modes with the choice s+ 
are orthogonal to those with s_ and so the two independent and properly normalized modes 
(positive frequency) are 

1 * '-^-^ <-*'>*•-"' (327) 



Consistency requires 



II C fc + = ; 1 r= i- k v) S - e~ ikv , 6a£ = , + * (-kr,) s - e~ ik \ (3.28) 
ca(r])yA.k a(ry)y4/c 



The behavior of modes in the coupled and decoupled cases are essentially the same when well 
within the horizon: compare the above to (1.2). We cannot modify the contribution to the 
spectral index from the subhorizon physics. 

For k <C aH the two modes must have the same power law dependence. Since in the 
krj — > limit (3.12) and (3.13) are a coupled system of Euler equations, we take = A' if 
and Sat = B'rf and substitute into (3.12) and (3.13) giving 4 , 

(6 - c 2 + p(c 2 - 2)){2B'h + A'(c 2 - 2)p) = 0, (3.29) 

and 

5'(24/i 2 - 4c 2 (/ 2 + h 2 )) + (c 2 - 2){-2A'c 2 h + B'p(6 - c 2 + p(c 2 - 2))) = 0. (3.30) 



4 We would like to thank Krzysztof Turzynski for pointing out a crucial oversight in an earlier version of 
this manuscript which modifies the previously stated conclusions. 
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From the first equation we have one solution with 



for which the second equation implies 



. _ A'c 2 {c 2 - 6)/i 

" ~2(-6h 2 + c 2 (f 2 + h 2 )Y (3 ' 32) 



Returning to (3.29), the remaining solutions have 

A'{c 2 - 2)p 



(3.33) 



and substituting into (3.30) gives 

p(c 4 (p - l)p + 4(6/i 2 + p(p - 3)) + 4c 2 (/ x 2 -f 2 -h 2 + 2p- p 2 )) = 0. (3.34) 
This has three solutions: one is p = 0, and the others are 



(c 2 -6)±V(c 2 -6) 2 -16c 2 A 

J* = . ( 3 - 35 ) 

where we have defined 

A = fl - f 2 - (1 - 6/c 2 )h 2 . (3.36) 

The case p = deserves special mention since in this case (3.29) and (3.30) imply that B' = 0. 
A more careful treatment shows that for this mode Cfc an d S^k do not behave as the same 
power for small rj and Ck = A'(l + 0(k 2 r] 2 )), Sat = A'(0 + 0(k 2 r] 2 )). This mode may be 
interpreted as the usual constant Ck mode that appears in single field models. Since Saj- 
vanishes faster that £fc, this mode is entirely adiabatic at long wavelengths. 

Allowing for all four powers and matching the fields and derivatives at horizon cross- 
ing we obtain the following approximate solution for the positive frequency modes (similar 
expressions apply for both modes I and II) 

C + = aiA; (( c2 - 6 )/( c2 - 2 )- 1 / 2 +P)(_ r/ )( C 2 -6)/( C 2 - 2 ) +a2fc (-l/ 2 +P) 

k ( - p -~ 1 / 2+ P\-ri) p - , (3.37) 

5 a + - aiC 2 (c 2 - 6)/i ((c 2 -6)/(c*-2)-l/2+p) ( _ x(c 2 -6)/(c 2 -2) /„ ogx 

k - 2(-6h 2 + c 2 (f 2 + h 2 )) k ( "> {S - S8} 



(c 2 - 2) 
2/i 



(0 + 0{k 2 r] 2 ))a 2 k^ 1 / 2+ p) + a 3 p + 1 / 2+ p^p+ + a 4 p_ 



where ai ~ a 2 ~ a 3 ~ 04 ~ ( — 7 ?o) p - As a result the exact adiabatic and isocurvature spectral 
indices of the modes that dominate at long-wavelengths are given by the minimum of 



! ^ ' ^2—^ > 4 " VV - 6) 2 - 16c 2 A J . (3.39) 
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The first entry corresponds to the constant £fc mode, the second to the solution described in 
(3.31) and (3.32), and the final term to the p± solution given in (3.33) and (3.35). 

It is now apparent that the third term in brackets can be made arbitrarily close to n s = 1 
by adjusting the field space metric and potential couplings to give the appropriate value of 
A. In particular exact scale invariance is achieved for 

A = - 1 -(c 2 -3) = ~w. (3.40) 

However, for this term to be dominant at long wavelengths we require the other two modes to 
give n s > 1. For the first term this implies c 2 > 2 or w > — 1/3. Consequently we immediately 
conclude that this last term can only be dominant and scale- invariant for contracting solutions. 
For the second term to be subdominant we require c 2 > 3 and so w > 0. 

In summary, there are three different ways of achieving scale-invariance, depending on 
which entry in (3.39) is dominant: 

• The first term in (3.39), at c 2 = 0, w = —1. This is an expanding, nearly de Sitter, 
inflationary universe. In this case the final term in (3.39) is also scale-invariant any A. 

• The second entry in (3.39), at c 2 = 3, w = provided also that A > 0. This is a 
contracting solution previously discussed elsewhere. 

• The final entry, for a contracting background with w > and A < 0. 

In the next section we shall see that the condition that the contracting scaling solutions are 
attractors is A > and w > 1. Therefore none of the contracting backgrounds that give 
scale-invariant adiabatic perturbation spectra are stable. 

The adiabatic power spectrum for the case where the third term dominates can be ex- 
pressed in a form similar to that for slow roll inflationary models 

™ - ?(^r is- <3 - 4i) 

up to a factor typically of O(l) 5 , where if* is the Hubble constant at horizon crossing, H en d 
the Hubble constant at the end of the generation process {i.e. reheating) and k p h ys is the 
physical wavenumber k/a(rjf) at the end of the generation process. The ratio of isocurvature 
to adiabatic power is given by 

P, ((^-6)- V(c2-6) 2 -16c2a) 2 

f = A w? -• {3A2) 

and so depends on the three parameters (c, A,/i). By taking fi to be suitably large whilst 
keeping (c, A) fixed we can put an arbitrary fraction of the power into the adiabatic spectrum 



3 It is feasible that this factor which depends on c, /i and A could be tuned to be arbitrarily large or small. 
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without affecting the adiabatic spectral index. Therefore it is not necessary to create large- 
amplitude isocurvature fluctuations in order to generate scale-invariant adiabatic fluctuations. 
In any event the amount of isocurvature modes that will be observed in the cosmic microwave 
background depends on the details of the reheating process, or its analogue in a contracting 
model. 

4. Stability 

We next determine whether the background scaling solutions are stable to small homogeneous 
perturbations. Solutions that are stable have the attractive property that they do not require 
their initial conditions to be fine-tuned. The stability properties are also intimately connected 
with the spectral indices of long-wavelength perturbations, since the stability analysis is 
cosmological perturbation theory in the k/a — > limit. The stability of selected two-field 
models, with some specific choices of field space metric and potential, has been studied before 



In our analysis the conditions required for scale-invariant long-wavelength perturbations 
will reappear as conditions on the stability of background solutions. These conditions can 
be conveniently expressed in terms of the parameters c and A that were introduced in the 
preceding section. In the expanding case, the scaling solutions that produce nearly scale- 
invariant adiabatic perturbations are always stable. In the contracting case, the relevant 
solutions are always unstable. 

To begin, we write the equations of motion for the homogeneous components of the fields 
as an autonomous dynamical system, and then investigate the stability of its fixed points. 
The original equations of motion are given in (3.2-3.4). For the time being we will assume 
Vq < and c > \/6. We then define 



These dimensionless variables, with the exception of z, express the fractional contribution 
of each term to the total energy density. The Friedmann equation (3.2) now appears as the 
constraint equation 



[27, 25]. 




(4.1) 




1. 



(4.2) 




(4.3b) 



(4.3a) 



dz 




(4.3c) 



din a 




(4.3d) 
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As required (4.3a-4.3d) propagate the constraint (4.2). There is a fixed point at e 
(xq, Vo, z , w ) = I —^= . 4 / — - 1 



hdf/da c 2 



67' V 6 / 



fdh/da 6/ 



,0,0 , 



(4.4) 



which corresponds to the scaling solution studied in previous sections. The values of xq and 
j/o express the requirement that a scaling solution must maintain a constant ratio between 
the potential and kinetic energy of the scalar field, which is itself determined by c. 

To determine whether the fixed point is stable, we expand around it, taking x = xq + 5x, 
etc. The constraint (4.2) implies that 



x Q - 



dx 



'dlna ^°'dlna : 

and so we solve the constraint and eliminate Sy as an independent variable by setting 



dy 



5y = 



5x. 



(4.5) 



(4.6) 



With this we have three equations for the perturbations, consistent with the constraint, which 
are 

c 2 - (i 



dx 
din a 
dw 

din a 

dz 

din a 



2 - fa "V^ + 24 fc ' 



6c/i 
-ox 



2 6 



Sw — r5z, 



= V6Sw, 



where we have defined 



36 



c\ft -h~ h 2 ) - ~^h 2 - 3(ft - 2[f 2 + 2h 2 ]) 



V6(c 2 - 6) 

The equations of motion for the perturbations are in the form 



d5xj 



Mj k Sx k , 



(4.7a) 
(4.7b) 
(4.7c) 

(4.8) 
(4.9) 



din a 

and so the stability of the fixed point is entirely determined by the eigenvalues of Mj h . These 
are 



An = 



c 2 -6 



c 2 - 6 



1± 4/1 



16c 2 



(c 2 - 6) s 



(4.10a) 
(4.10b) 



6 Note if / jCT = there are additional fixed points at (±1,0, zo,0) with z constant, however these are 
physically uninteresting. 
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where as before 

A = f! - h 
The eigenvectors corresponding to A± are 



-2cf ± 5x + (c 2 - 6) 



1± 



16c 2 



(4.11) 



(4.12) 



y (c 2 - 6) 2 

and the eigenvector corresponding to Ao is 

[-3{-6h 2 + c 2 [f 2 + h 2 \) - c 2 {c 2 - 3)A] Sx + 3(c 2 - 6)5w + 6>/6c/i<fc. (4.13) 

For a contracting solution, in which a — ► 0, all of the eigenvalues should be positive to ensure 
stability. The expression for Ao thus requires 

c > \/6 stable, contracting case, (4-14) 

in accord with previous results with a single scalar field [5]. With this, A + is always positive, 
but A_ is positive only when 



A > stable, contracting case. 



(4.15) 



To study the expanding case, one must start from the original equations of motion but assume 
Vq > 0, as required by the scaling solutions with c < \/6- The corresponding autonomous 
system, and constraint, can be derived more simply from the analytic continuation y — > iy in 
(4.2) and (4.3a-4.3d). The eigenvalues of the matrix Mj k are precisely the same, but since 
a is increasing we require the eigenvalues to be negative in order to obtain a stable solution. 
This results in the conditions 

c < \/6 stable, expanding case, (4.16) 

again in accordance with results for a single scalar field. In addition the expression for A_ 
requires that 

A > stable, expanding case, (4-17) 

to ensure stability of the expanding solution. 

Therefore, by comparing these results to those of the previous section, we see that the 
expanding universes that produce nearly scale-invariant adiabatic spectra are always stable 
background solutions. Furthermore, the contracting universes are always unstable. Within 
this framework it is also possible to study contracting universes with f\ = that produce 
scale-invariant isocurvature perturbations that are later converted into adiabatic modes. One 
must carefully redo the analysis since, as described in the last section, the decoupled /i = 
case is not necessarily the limit f\ — ► of the coupled case. In the end, one finds that 
the stability properties are controlled by the analogue of A with f\ set to zero. Since the 
isocurvature spectral index is smooth in the f\ — ► limit, by comparing the results of this 
section to those of the previous one, we can conclude that contracting background solutions 
that produce a scale-invariant spectrum in isocurvature modes alone are unstable. Whether 
this is a problem depends on the specific model one considers. 
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5. Conclusions 



In the contracting case, we have found that a nearly scale-invariant spectrum of perturbations 
in the variable £ can be generated continuously as modes exit the horizon for a one-parameter 
family of (c, A). The problem of generating such a spectrum in models where primordial per- 
turbations are laid down during a contracting phase has posed a significant challenge, as 
well as sparking lively debate. Other proposed mechanisms to generate the proper spectrum 
involve a two-stage process, by which isocurvature perturbations with the proper spectral 
index are generated on horizon exit, and then later converted into adiabatic perturbations. 
The mechanism we describe here is more akin to the classic inflationary mechanism in which 
the desired density perturbation is created in a single step as modes exit the horizon. Fur- 
thermore, this can be achieved entirely within a four-dimensional effective theory, without 
recourse to higher-dimensional effects or higher-derivative interactions. There remains the 
challenge of a smooth big crunch/big bang transition, but assuming that this is possible, it 
seems that within the context of a simple model it is certainly possible to generate the re- 
quired spectrum of primordial density perturbations in cosmological models with a collapsing 
phase. 

In the expanding case, we find the expected conclusion that a nearly scale-invariant 
perturbation spectrum can only be obtained for w « — 1, i.e. close to slow roll, regardless of 
how we tune the other couplings (/i, A). This is because there always exists a constant mode 
in Q which necessarily dominates the long-wavelength power spectrum and whose spectral 
index is solely determined by the equation of state w. 

In this work we have obtained these results by studying idealized two scalar systems that 
admit scaling solutions. The scaling symmetry is powerful enough to enable the background 
evolution and the spectral indices of perturbations to be computed exactly. The scaling sym- 
metry is also general enough to provide reasonable representatives of a large class of models 
in which the equation of state parameter w varies slowly during the time that primordial per- 
turbations are laid down. Within the context of these scaling solutions it has been possible 
to completely characterize the conditions under which a scale-invariant spectrum of density 
perturbation is created, after the spirit of previous studies in the one-field case [67, 5, 68]. 

There are several directions in which the work presented here should be extended. With 
respect to model building, the possibility of achieving a scale-invariant adiabatic perturbation 
spectrum with a simple two-derivative action opens up new avenues for realizing cosmological 
models with a contracting phase. This is especially promising with respect to string theory, 
where multiple scalar fields and curved field spaces are common. Thanks to the instability 
in the contracting solutions, complete cosmological models that incorporate this mechanism 
must explain how the scalar fields begin with the proper initial conditions in field space. 
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